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Abstract. A simple version of the q-deformed calculus is used to generate a pair of 
q-nonlocal, second-order difference operators by means of deformed counterparts of Dar- 
boux intertwining operators for zero factorization energy. These deformed non-local 
operators may be considered as supersymmetric partners and their structure contains 
contributions originating in both the Hermite operator and the quantum harmonic oscil- 
lator operator. There are also extra ±x contributions. The undeformed limit, in which 
all q-nonlocalities wash out, corresponds to the usual supersymmetric pair of quantum 
mechanical harmonic oscillator Hamiltonians. The more general case of negative factor- 
ization energy is briefly discussed as well. 
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I present herein a simple g-deformed procedure |jj for the basic case of the one- 
dimensional quantum harmonic oscillator, by which I build a 'supersymmetric' pair 
of q nonlocal operators possessing terms whose q — ► 1 limits belong to either Hermite 
polynomial operator or Schroedinger quantum oscillator operator. There are ±x 
extra terms as well. The procedure is based on the idea of using, as fundamental 
tools, a sort of deformed counterparts of the intertwining operators encountered in 
the area of Darboux transformations 0. I shall use their factorization property 
to get the q-deformed second-order operators which, being q non-local, may be 
considered as more general than both the usual Hermite one and the quantum 
mechanical harmonic oscillator operator. 

The standard Hermite operator Oh reads (D = d/dx) 

6 H = D 2 - 2xD + 2n , (1) 

and gives rise to the equation for the Hermite polynomials H n (x), OjiH n (x) = 0. 
Writing 6 H = D 2 -2xD-2+2(n+l), I shall treat D 2 -2xD-2 as the Fokker-Planck 
(FP) part of the Hermite operator for a stationary transition-probability density, 
since —2xD corresponds to the (^r)D drift contribution in the FP stationary op- 
erator (drift potential U = —x 2 ), whereas —2 stands for the d 2 U/dx 2 contribution 
of the FP drift. The last term 2(n+ 1) gives the departure of the Hermite operator 
from the corresponding FP stationary operator for which polynomial oscillations 
are not allowed, and in fact is responsible for turning the FP interpretation into a 
formal one and not a physical one. As well known for this basic case, by means of 
the functions <f> n — e~ x / 2 H n (x) one can go to the operator 

O = -D 2 + [x 2 - (2n + 1)] , (2) 

which, in the <p n space, is essentially the Schroedinger quantum harmonic oscilla- 
tor operator up to a scaling, choose-of-units factor. One should notice that this 
usage of the <fi n functions leads to the loss of one half of the 45 drift contribu- 
tion. The remaining half gets the famous zero-point energy interpretation when the 
scaling \0 < p is performed. In the FP interpretation, the latter scaling corresponds 
to setting the diffusion constant equal to 1/2 and provides the usual quantum me- 
chanical harmonic oscillator wavefunctions N n <j) n , where N n = (2™n!v / 7r) -1 / 2 is the 
normalization factor. 

I now briefly recall that in the case of the one-dimensional Schroedinger oper- 
ator within the context of supersymmetric quantum mechanics (SUSYQM) || the 
standard Darboux transformation operator reads 

T = —t u {x) + D = —u'(x)/u(x) + D , (3) 

where the prime denotes the derivative with respect to x. When acting on the solu- 
tions tp n {x) of the initial Schroedinger equation hotp n (x) = E n tfi n (x), it transforms 
them into the solutions of another Schroedinger equation hi<p n (x) — E n ip n (x), 
fn(x) = N n Tip n (x), with the same eigenvalues E n . Henceforth, I will put the 
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ground state energy equal to zero, Eq — 0, since this does not affect in any way 
the results. The new exactly solvable Hamiltonian has the form hi = ho + AV(x), 
where the potential difference is of Darboux type AV(x) = — 2(ln it)". The func- 
tion u = u(x) is a so-called transformation function, being a solution of the initial 
Schroedinger equation hou(x) = eu{x), with e < usually known as the factoriza- 
tion energy. It is well established that when e < one can work with a nodeless 
transformation function by performing an analytic continuation Thus, u(x) ^ 
for any value of the variable and \/u{x) is not a square intcgrable function. In this 
case u ^ Hi and the set {| (p n )} is a complete basis in the Hilbert space Hi provided 
the initial system {| ip n }} is complete. The operator T + = —t u (x) — D provides 
the backward transformation \ip n ) = N n T + \(p n ), and together with T allows for the 
following factorizations 

T+T = h -e, TT + = hi - e . (4) 

The operators T and T + are well defined VV> S Hi and are conjugated to each other 
with respect to the inner product in the Hi space. 

My purpose now is to get g-deformed second-order operators by means of de- 
formed counterparts of the aforementioned intertwining operators. I still have to 
present some definitions and rules of the deformed calculus. Since the independent 
variable is maintained commutative, the employed version of the deformed calculus 
is similar to that previously used by some authors to deform the Coulomb problem 
||. Symmetric definitions of the q- number [x] q — q Z q -i an d ^-derivative 



f( qx) - fjq^ x) 



dm = jyq :r j ^ ' (5) 



are used together with some basic rules of Jackson's calculus [Q such as D q 



-l r>2 



DgX n = [n} q [n - l] q x n - 2 , D q (FG) = (D q F)G(qx) + F^x^D^) for 
any two functions F and G, respectively. The definition of the ^-exponential is 

~ x n 

^)-Y.^ r (6) 

n=0 1 lq 



which reduces to the usual exponential function as q — > 1, and moreover is invariant 
under q — > q^ 1 . 

The main idea of this work is based on the following scheme. First, to employ 
as Darboux transformation functions deformed counterparts of the oscillator vacua 
i/jq oc e q (/3x 2 ), where (3 — ±1/2 for the irregular and regular vacuum, respectively. 
Second, to exploit the factorization property of first-order deformed operators of 
the form 

Tl = D q -^L = D q -/3 q (x 2 )x , (7) 

Vq 

Tl = -D q -^ = -D q -l3 q (x 2 )x, (8) 

Vq 
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where 

^ = p ( geg ( g ^) + g -^(g-^) \ (g) 
V e q (px ) J 

The form of f3 q {x 2 ) is a result of Jackson's calculus rules. As one can see, the 
T? and T£ operators have been written by analogy to the continuous intertwining 
operators. A straightforward calculation gives the second-order deformed operators 
that can be obtained from the products T^T? and TlT£, respectively. One gets 



01 = TlTl = -^-[(AftJ^l+^^J^l^+bp^^x^il^+Iftl.;- 2 ! 2 )]^,, 
and 



' 2 ™ 2 ) 
(10) 



6) = TlTl = -^ 2 +[(A/3 9 )^] + [/3 9 2 (x 2 ) a; 2 ]^-[ g (i) g /J 9 ( 2 : 2 ))a ; ]^-[/3 9 ((7- 2 a ; 2 )]^ ;l 

(11) 

where 

= Pq{x 2 ) - q~ 1 {3 q (q~ 2 x 2 ) . (12) 

The operators 0\ and may be considered as supersymmetric partners since they 
have been built according to the well-known SUSYQM method. At this point one 
should notice the interesting mixed structure of the two q non-local operators that 
entail parts of both Oh and O^. The directional subindices indicate the argument of 
the solution on which the nonoperatorial parts act. The two operators are nonlocal 
operators whose space of solutions are the functions <j)n\x) oc e q {—x 2 /2)Hn\x), 
where the deformed Hermite polynomials can be defined by a q deformed Rodrigues 
representation 

ffW = (-l) n e q {x 2 )D^(e q {-x 2 )) . (13) 

Notice that the D q terms in (10) and (11) are identical but opposite in sign and 
correspond to the first derivative drift term in the Hermite differential operator. Of 
course, if one prefers the FP interpretation the two operators should be multiplied 
by (-1). Writing the finite difference x(q — q^ 1 ) = xAq = A q x, which for q — > 1 
is assumed to be a q scaling way of going to the infinitesimal limit dx, the q drift 
parts go to zero, whereas in the same limit the potential and zero point sectors take 
forms identical to those of the undeformed case. More precisely, the undeformed 
limits read 

61 = h Q = -D 2 + Plx 2 + fa (14) 

and 

6) = hi = -D 2 + 0lx 2 - (3i . (15) 

(14) and (15) are the usual quantum mechanical supersymmetric partner Hamilto- 
nians for this case. 

In SUSYQM terminology, only the case of zero factorization energy e = has 
been tackled up to now, but following a suggestion of Bagrov and Samsonov 
there is no difficulty to sketch the procedure for the more general case e < 0. First, 
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the deformed Schroedinger solution corresponding to the excited harmonic oscillator 
states can be sought in the form 

^\x)<xH^{x)e q (-x 2 /2) . (16) 

Next, in order to avoid any singularities, it is convenient to perform an i-rotation 
x — > ix, leading to 

u^(x)(xH^(ix)e g (x 2 /2), p = 0,1,2,3... (17) 

The undeformed functions are solutions of hoiip — — (p + l)u^ and are 
nodeless on the full line for even p = 2k. Therefore, they have been used by Bagrov 
and Samsonov as Darboux transformation functions to generate a family of regular 
potentials, which, according to an interpretation due to Veselov and Shabat , has 
a spectrum made up of 2k + 1 segments with equidistant levels. This immediately 
suggests using u^\x) for even p = 2k as Darboux transformation functions in 
the deformed case. Thus, the intertwining operators can be calculated according 

to T| fc , = D q q - ! ^ k - and T$ k _ = —D q i^yS an d again by exploiting the 

U 2k _ ' U 2k 

factorization property one is led to second-order deformed, non-local operators of 
more complicated formulas than (10) and (11) for which they are not written down 
here. 

In conclusion, a pair of q non-local second-order g-differential (g-difference) op- 
erators have been introduced in this work by means of a particular g-deformed 
intertwining based on g-deformed oscillator vacua as Darboux transformation func- 
tions. These operators display a mixed structure between the Hermite operator, 
to which they are similar as regards the first derivative term, and the quantum 
mechanical oscillator operator, to which they are similar as regards the x 2 potential 
and zero-energy contributions. On the other hand, they present a supplementary 
q non-local ±x potential contribution with no counterpart in either Hermite poly- 
nomial operator or Schroedinger x 2 oscillator operator. All these features suggest 
many possible applications, e.g., in mesoscopic physics. A more general case corre- 
sponding to negative factorization energies of the type e m = — (to + 1), where to is 
an even positive integer, has also been briefly described. 
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